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INTRODUCTION 


The  problem  concerning  a  point  explosion  in  a  gas  consists,  as 
is  known  [1,  2],  in  determining  the  parameters  of  the  gas  flow,  which 
has  been  caused  by  the  Instantaneous  release  of  finite  energy  E°  at 
a  point  (spherical  explosion),  along  a  straight  line  (cylindrical 
explosion)  or  along  a  plane  (flat  explosion).  In  the  cylindrical 
case  tl’.e  energy  of  the  explosion  E°  is  calculated  per  unit  length, 
and  in  the  flat  case  —  per  unit  area. 

The  problem  concerning  a  point  explosion  is  one  of  the  basic 
transitional  problems  of  the  mechanics  of  a  continuous  medium.  It 
can  be  formulated  and  examined  both  for  various  ideal  compressible 
mediums  and  also  for  continuous  media,  in  which  It  is  impossible  to 
disregard  viscous  stresses.  When  studying  this  problem  for  a  gas, 
in  a  number  of  cases  one  should  take  into  account  influence  of  various 
physical  and  chemical  parameters  and  phenomena  on  the  flow.  Thus, 
it  is  possible  to  take  into  consideration  the  variability  of  the  initial 
density  of  the  gas,  to  investigate  the  effect  of  electromagnetic 
fields  on  the  motion  of  an  electroconductive  gas,  to  calculate  the 
chemical  reactions  during  an  explosion  in  a  combustible  mixture,  and 
to  examine  the  processes  of  equilibrium  and  nonequilibrium  dissociation 
and  ionization.  The  questions  mentioned  above  in  the  point  explosion 
theory  are  now  being  developed  by  various  authors . 
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This  article  is  dedicated  to  the  problem  of  a  point  explosion 
with  counterpressure  in  a  quiescent  nonviscous  and  nonheat-conducting 
ideal  gas  with  the  constant  adiabatic  exponent  y.  Pressure  p,,,  and 
density  p#  of  the  undisturbed  gas  are  considered  constants.  As  a 
result  of  numerical  solution  of  this  problem,  tables  of  the  gas- 
dynamic  functions  of  a  point  explosion  are  calculated. 

This  article  presents  tables  for  the  flat,  cylindrical,  and 
spherical  explosions  in  an  ideal  gas  for  the  values  of  y  equal  to 
1.3,  1.4,  and  5/3  (for  y  *  1.3  the  data  are  given  only  in  the 
cylindrical  and  spherical  cases).  The  tables  contain  the  basic 
functions  describing  the  field  of  the  gas  flow  over  a  wide  time  span. 
Also,  some  possible  applications  of  these  tables  for  solving  other 
problems  in  the  gas-dynamics  of  explosion  and  in  hypersonic 
aerodynamics  are  explained  in  the  introductory  section. 

The  literature,  already  contains  tables  for  an  explosion  in  a 
gas  with  counterpressure  [3]-  However,  these  tables  are  only  for 
the  spherical  case  when  y  =  1.4.  Let  us  note  also  that  we  [4]  first 
calculated  the  tables  of  the  gas-dynamic  functions  for  the  initial 
stage  of  a  point  explosion,  i.e.,  for  sufficiently  large  pressure 
gradients  at  the  shock  wave  front.  The  tables  now  being  published 
are  the  natural  continuation  of  the  tables  for  the  initial  stage  of 
an  explosion. 

The  calculation  of  the  current  tables  of  gas-dynamic  functions 
was  accomplished  on  BESM-2  and  BESM-3M  computors  of  the  Computation 
Center  of  the  USSR  Academy  of  Sciences.  The  numerical  method  of 
solving  explosion  problems,  used  tb  calculate  the  tabulated  data, 
is  comprehensively  described  in  works  [5,  6],  They  give  basic 
equations  and  a  number  of  useful  formulas  and  reductions,  which  are 
not  reproduced  here. 

Taking  the  opportunity,  the  authors  wish  to  express  their 
gratitude  to  Ye.  Bishimov  for  his  assistance  in  compiling  programs 
for  the  BESM-3M  computor ,  to  L.  S.  Bark,  who  supervised  the  checking 
of  the  tables  for  diversities,  and  to  V.  P.  Karlikov  for  his  valuable 
remarks  under  the  editing  of  the  manuscript. 
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1.  METHOD  OP  CALCUATING  DIMENSIONLESS 
GAS -DYNAMIC  FUNCTIONS. 

The  functions  of  the  explosion  problem  to  be  determined  (pressure 
p,  density  p,  velocity  v,  temperature  T,  total  energy  per  unit 
volume  E)  depend  upon  geometric  coordinate  r  and  time  t.  Coordinate 
r  in  the  plane  case  (v  *  1)  represents  the  distance  from  the  plane 
of  the  explosion,  in  the  cylindrical  case  (v  ■  2)  —  the  distance  from 
the  axis  of  tne  explosion  and  in  the  spherical  case  (v  *  3)  -  the 
distance  from  the  point  of  the  explosion.  The  law  of  motion  for  a 
shock  wave,  i.e.,  the  dependence  of  the  shock  wave  coordinate  rR  upon 
time  t,  will  also  be  determined  when  the  problem  is  solved.  The 
functions  specified  above  depend  parametrically  upon  the  adiabatic 
exponent  y,  the  energy  of  the  explosion  E°  and  the  initial  values 
of  the  pressure  pM  and  the  density  pw. 

Various  numerical  methods  are  used  for  solving  the  problem  of  a 
point  explosion.  Thus,  in  the  above-mentioned  work  [3],  a  special 
numerical  method  based  on  a  scheme  of  finite  differences  was 
developed.  The  calculation  of  the  gas-dynamic  functions  contained  in 
the  current  tables  was  done  by  the  integral  relationship  method. 

The  problem  is  easily  solved  in  dimensionless  variables.  Let  us 
introduce  the  following  independent  variables: 
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where  c  -  the  velocity  of  the  shock  wave;  a,,,  -  the  speed  of  sound  in 
a  quiescent,  which  is  expressed  thus: 


9 

fTf 


The  range  of  change  for  the  variables  E  and  q  is  giveO  by  the 
Inequalities  0  <  ?  <  1,  0  <  q  <  1.  In  this  range  we  will  single  out 
the  central  interval  limited  by  lines  E  ■  0  and  E  ■  EQ(q)>  in  which 
the  solution  is  found  with  the  help  of  the  asymptotic  formulas 
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♦  0\* 


1»»  -  i\ 
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where 
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v  ♦  i(r-» 
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Here  p0»  p(t0*  >).  Po*  •>.  •  »(t0,  f)  represent  the  values  of  the 

pressure,  the  density  and  the  velocity  at  the  boundary  of  the  central 
interval,  i.e.,  on  linO  E  ■  Eq (q ) *  To  select  a  vplue  for  Eq»  the  i 

Lagrangian  coordinate  of  the  particle  is  fixed  so  that  ih  the  initial 

‘  ,  ; 

stage  of  the  explosion  it  is  possible  to  determine  the  pressure  with 
acceptable  acouracy  according  to  the  appropriate  asymtotic  formula 
from  ( 2 ) .  ' 


The  solution  in  the  domain  located  between  the  lihe  E  *  Eq(9) 
and  the  shock  wave  E  *  En  ■  1  is  formed  by  the  method  of  integral 
relationships.  In  the  n-th  approximation  this  domain  (Fig.  1)  is 
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divided  into  n  strip?  by  n  -  1  the  intermediate  lines: 


i  ^  i (9) m  ^  ^  ^  li  (*•  1.2.  1).  (3) 

,  I 

For  convenience  the  new  dimensionless  unknown  functions  are  introduced. 
The  Initial  system  of  gas-dynamics  equations  is  converted  to  dimension¬ 
less  variables  and  Is  reduced  to  divergent  form.  Every  equation  of 
the  obtained  system  is  integrated  with  respect  to  £  from  the  value 
5  -  (Z  «  0,  1,  ....  n  -  1)  to  the  shock  wave  5  -  £n  -  1,  which 
leads  to  a  system  integral  relationships .  The  approximation  of  the 
intergrands  by  Lagrange  interpolation  polynomials  with  interpolation 
points  on  lines  £  *=  £fe(q)  (fe  -  0,  1,  ...»  n)  gives  an  approximating 
system  of  ordinary  differential  equations  with  respect  to  q  for 
determining  the  unknown  gas-dynamic  functions  on  lines  £  ■  £j(q)  and 

the  dimensionless  coordinate  Rn  of  the  shock  wave. 


The  initial  equations  and  conditions  of  problem  also  make  it 

possible  to  derive  differential  equations  for  determining  the  value 

of  and  of  dimensionless  time  x  as  functions  of  the  variable  q 

Introduced  above .  Time  t  and  the  shock  wave  coordinate  rn  are 

connected  with  the  values  x  and  R  in  the  following  manner: 

n 
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where  dynamic  length  r°  and  time  t°  are  expressed  thusly: 


For  a  complete  approximating  system  of  differential  equations 
the  Cauchy  problem  is  numerically  solved  with  the  initial  data  at  a 
certain  small  value  q  ■  qQ  taken  from  the  tables  in  [4],  which  contain 
the  solutions  of  self-simulating  and  linearized  explosion  problems. 

The  tables  presented  here  were  calculated  by  dividing  the 
lntergration  range  into  the  eight  strips,  i.e.,  when  n  ■  8,  while 
in  the  case  of  constructing  an  approximating  system  of  ordinary 
differential  equations,  the  integrads  in  the  integral  realizations 
of  the  problem  were  replaced  by  two  4th  order  joining  interpolation 
polynomials.  Besides  the  basic  calculation,  when  n  *  8,  calculations 
of  lower  approximations  when  n  ■  1,  2  and  4  were  also  made;  in  all 
these  cases  different  forms  of  approximating  polynomials  were  used. 
Comparison  of  the  obtained  results  characterized  the  accuracy  and 
practical  convergence  of  the  method  of  integral  relationships  in 
the  explosion  problem. 

To  check  the  accuracy  of  the  numerical  solution,  the  integral 
laws  of  conservation  of  mass  and  enprgy  in  the  total  volume  of  the 
moving  gas  were  also  used.  Let  us  designate  by  the  relative 
mass  error,  which  is  equal  to  ratio  of  the  difference  between  the 
calculated  mass  of  the  moving  gas  and  the  initial  mass  of  gas  in  the 
volume  limited  by  the  shock  wave  front  to  the  initial  mass  of  the  gas 
in  this  same  volume.  We  will  designate  as  e2n  the  ratio  of  the 
difference  between  the  calculated  total  energy  of  the  perturbed  gas 
and  the  initial  total  energy  in  the  previously  mentioned  volume  of  gas 
to  the  Initial  total  energy  of  the  gas  in  this  same  volume.  The 
quantity  e2n  gives  the  relative  energy  error. 


Furthermore,  for  initial  value  q  ■  qg,  the  error  was  checked 
by  the  function  ip  -  (qp/p,,,)  1/y,  conditioned  by  the  selection  of  the 
central  interval  £g(qg).  The  quantity  e^n  is  introduced  in  the 
following  manner: 

.  WW-tIM.) 

**  y(0.90) 

During  the  assignment  of  the  initial  data  the  values  qQ  and 
£g(qg)  were  selected  so  that  the  conditions 


<  6  j«  • 

would  be  satisfied  at  small  q.  Here,  e°n,  e°n,  e°n  -  assigned  numbers, 
whose  selection  is  determined  both  by  the  requirement  for  their 
smallness  and  by  the  conditions  of  the  proximity  of  the  numerical 
solution  to  the  linearized  solution  in  the  vicinity  of  q  ■  q^  and 
the  absence  in  the  region  of  the  small  q  of  a  physically  unreal 
oscillation  in  the  solution  at  large  n  (n  >  4). 

In  the  conducted  calculations ,  the  value  of  qQ  was  always  taken 

as  equal  to  0.05-  The  values  of  the  initial  errors  e.  ,e_,e_ 

In*  2n*  3n 

for  all  the  calculated  cases  are  shown  in  the  following  table. 


V  1 

•  1 

v  *  2. 

v-  3 

Y 

1.4 

S/S 

1.3 

1.4 

5/3 

1.3 

■sa 

Fflm 

■in 

0,043 

0,007 

0,014 

0*009 

0,003 

0,016 

0,009 

0,021 

-0,002 

-0,001 

-0,002 

0,006 

-0,002 

0,003 

*3* 

0,078 

0,028 

0,023 

0,025 

0,013 

0,032 

0,026 

0,078 
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Let  us  note  that  the  errors  in  the  initial  da^a  fade  rather 
rapidly  during  calculation.  . 


Internal  checking  of  the  accuracy  of  the  calculation  scheme  was 
additionally  accomplished  by  computing  the  dimensionless  pressure 
using  the  equation  of  conservation  of  total  energy  and  by  using  the 
equation  of  the  conservation  of  entropy  in  a  particle  of  gas.  The 
relative  divergence  in  the  pressure  values  determined  by  two  such 
methods  was  less  than  0.01. 

Curves  of  as  functions  of  q  for  n  ■  8  at  various  y  and  with 
v  *  1  (solid  line),  v  *  2  (crosses),  and  v  ■  3  (circles)  are  presented 
in  Fig.  2.  As  is  evident,  the  Integral  error  in  the  basic 
component  of  the  interval  of  change  comprises  from  0.002  to  0.02.  As 
for  the  absolute  value  of  e2n*  ^  was  max^-mum  at  the  beginning  of 
the  calculation,  i.e.,  when  q  ■  qQ,  and  then  it  dropped  to  values 
close  to  zero.  One  should  still  keep  in  mind  that  the  major  character¬ 
istics  of  an  explosion  (the  law  of  motion  for  the  shock  wave  and  the 
pressure  distribution  near  the  center  of  the  explosion)  are  derived 
with  greater  accuracy  than  the  field  of  the  gas  dynamic  functions.  In 
the  region  of  extreme  values  of  q  the  accuracy  of  the  calculated 
solution  is  reduced.  However,  for  the  initial  stage  of  the  explosion 
(q  <_  0.10)  it  is  possible  to  use  the  linearized  solution  in  [4],  and 
with  long  times  (q  >  0.90)  the  asymptotic  laws  of  the  fading  of  shock 
waves  can  be  used. 


The  evaluations  of  accuracy  indicated  above  made  it  possible 
in  the  greater  part  of  the  tables  to  present  the  values  of  the  functions 
to  three  decimal  places;  in  certain  cases  the  third  place  is  not 
completely  accurate. 
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2.  THE  STRUCTURE  QP  TABLES  AND  FORMULAS 
FOR  DETERMINING  THE  PHYSICAL  PARAMETERS 
OF  AN  EXPLOSION 

The  tables  contain  the  numerical  solution  of  the  problem  of  a 
point  explosion  in  an  ideal  gas  with  counterpressure.  They  encompass 
the  plane  (v  *  1),  cylindrical  (v  *  2)  and  the  spherical  (v  ■  3) 
cases  at  adiabatic  exponent  values  of  y  *  1.3  (only  for  cases  v  ■  2 
and  v  =  3),  y  =  1.4  (air)  and  y  ■  5/3  (monatomic  gas).  The  values 
of  the  basic  gas-dynamic  functions  are  given  as  functions  of  the 
variable  q  which  characterizes  time. 
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The  argument  of.  q  varies  from  q  ■  0.05  to  q  ■  0.90  with  constant 
step  of  Aq  -  0.05.  With  q  ■  0.05  and  q  •  0.10  data  sire  given  for  the 
linearized  solution.  The  values  when  q  •  0  can  be  obtained  from  the 
tables  in  [4]  for  the  self-similar  solution. 

At  the  beginning  of  the  tables  dimensionless  time  t,  the 
dimensionless  shock  wave  coordinate  Rn  and  the  boundary  of  the  central 
interval  depending  upon  q  are  presented.  The  following  relative 
functions,  selected  from  considerations  of  the  scale  of  the  tabulated 
data,  are  given  in  the  same  place: 


These  functions  give  the  values  for  pressure,  density,  velocity, 
temperature  and  energy  directly  behind  the  shock  wave  front.  The 
latter  values  do  not  depend  upon  v  and  are  expressed  by  these  formulas 


P  .UL--U.-.  Hp  .  p  .  y+  »  , 

"  <Y  ♦  I)*  P"  Y-  1  ♦  2q  P**  * 


v  «-2£-(|_ii)-  t  -  P"  •  f 

*  Yt  1  ’  •  »  RrP(|  '  c 


Y-l 


P  V 

«_  +  *"» 


(6) 


In  accordance  with  definition  (1),  the  velocity  c  of  the  shock  wave 
will  be 


In  the  formulas  in  (6)  Rp  designates  the  gas  constant. 

The  basic  contents  of  tables  are  a  description  of  the  fields  of 
the  gas-dynamic  functions ,  which  for  the  sake  of  convenience  are 
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related  to  the  corresponding  values  of  (6)  and  (7)  at  the  shock  wave 
front.  For  the  pressure,  density,  velocity,  temperature  and  energy 
the  following  relative  variables  are  Introduced: 


The  tables  give  the  values  of  these  functions  on  lines  £  ■  £^(q)  (7  * 

0,  1,  ...,  n  -  1),  i.e.,  at  the  boundary  of  the  central  interval 
£  *  £Q(q)  and  on  the  intermediate  lines  £  *  ^(q)  determined  by  equation 
(3).  The  values  of  the  functions  at  points  £  that  do  not  coincide 
with  the  end-points  £  *  £^  are  found  by  interpolation.  Here  it  is 
advantageous  to  use  an  interpolation  polynomial  not  of  the  power 
n  ■  8,  but  of  smaller  powers  (for  instance,  the  second  or  third), 
while  conducting  local  Interpolation  at  the  nearest  nodal  points, 
which  is  entirely  sufficient  in  terms  of  the  accuracy  of  the  calculated 
solution. 

Let  us  note  certain  peculiarities  in  the  calculation  of  gas- 
dynamic  functions  inside  the  central  interval  when  0  <  £  <  £Q. 

Since  at  the  center  of  the  explosion  when  £  ■  0,  the  velocity 
and  the  density  are  equal  to  zero,  to  compute  the  values  v/c  and 
p/Pn  on  the  segment  [0,  £p]  it  is  possible  to  either  use  quadratic 
interpolation  over  the  known  values  of  these  functions  at  the  points 
£  *  0,  £  «  £q  and  £  ■  £1  or  use  the  asymptotic  formulas  in  (2). 

As  for  the  pressure,  its  precise  value  at  the  center  of  the 
explosion  is  unknown,  and  here  several  methods  of  determining  this 
value  approximately  are  possible.  Let  us  introduce  the  designation 
h  ■  p/pn,  where  this  ratio  of  the  pressures  when  £  *  0  we  will  note 
as  h«  *  h(0,  q). 

For  rough  estimation  in  conformity  with  (2),  it  is  possible  to 
assume  that  hg  ■  hQ.  To  find  the  pressure  inside  the  central  interval 
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it  is  also  possible  .to  use  standard  extrapolation  over  the  values  of 
h  at  the  points  5  ■  5Q  and  5  ■  However,  with  a  large  value  of 

the  central  interval  this  method  gives  the  value  of  the  pressure  at 
5*0  with  a  high  degree  of  error.  A  more  precise  extrapolation  for 
calculating  h  in  the  center  of  the  explosion  cam  be  conducted  in 
oonformity  with  asymptotics  in  (2)  using  the  following  formula: 


where  hQ  1  is  the  value  of  h  at  the  point  5  ■  belonging  to 

the  segment  [£q,  5-^]. 


The  determination  of  dimensional  physical  parameters  with  assigned 
specific  values  of  the  explosion  energy  E°  ,  the  pressure  p,,,  and  the 
density  p^,  of  the  undisturbed  gas  is  done  from  the  tabulated  data  without 
without  special  labor.  Time  t  and  the  shock  wave  coordinate  rn  are 
found  from  formulas  (4)  and  (5),  while  the  geometric  coordinate  r 
of  any  point  is  expressed  through  5  in  the  following  manner: 


r 


The  basic  physical  parameters  in  dimensional  form  are  calculated 
according  to  the  simple  relationships 
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where  the  functions  presented  in  the  tables  are  shown  in  brackets. 


It  is  necessary  to  further  note  that  the  numerical  solution  of 
the  explosion  problem  changes  quite  smoothly  as  a  function  of.  value  of 
the  adiabatic  exponent  y  in  the  range  of  the  examined  values  of  y. 

This  fact  makes  it  possible  with  the  help  of  interpolations  to  use 
the  given  tables  for  other  cases,  in  which  the  adiabatic  exponent 
y  is  different  from  the  tabulated  values. 

3.  CALCULATION  OP  CERTAIN  ADDITIONAL 
CHARACTERISTICS  OF  A  GAS  PLOW 


The  given  tables  can  also  be  used  for  computing  certain  supple¬ 
mentary  parameters  of  a  gas  flow  during  an  explosion,  which  are  not 
directly  incorporated  in  these  tables.  Let  us  present  a  number  of 
examples  of  such  an  application  of  the  tables. 

A.  Determination  of  Flow  Parameters  at 
Fixed  Points  of  Space 

Calculating  the  change  in  flow  parameters  at  fixed  points  in 
space  is  done  in  the  following  manner.  Suppose  that  a  point  is 
examined  that  has  a  fixed  dimensionless  geometric  coordinate  R  *  R#, 
where  R  *  r/r° .  The  perturbed  motion  of  the  gas  at  this  point  is 
initiated  from  that  instant  (or  the  corresponding  values  q„),  when 
the  shock  wave  arrives  at  this  point,  i.e.,  when  R  (q.)  *  R*.  For 
all  subsequent  q  the  following  values  are  determined: 


For  these  using  the  tabulated  values  of  the  functions  at  the  end¬ 
points  £  *  the  unknown  functions  at  the  considered  point  with 

the  coordinate  R„  are  found  with  the  help  of  interpolation. 
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I 


B.  Determination  of  the  Pulses  of  Excess 
Pressure  and  Velocity  Pressure 

Having  assigned  a  certain  value  for  coordinate  Rt,  using  the 
procedure  described  above,  let  us  find  in  this  point  of  space  the 
value  of  the  dimensionless  pressure  P  ■  p/p^  for  a  sufficiently 
large  number  of  instants  t.  In  terms  of  the  found  values  of  P(R§,  t)  * 
■  P„  It  is  possible  to  compute  the  complete  dimensionless  pulse  of 
the  excess  pressures 


S  (P.-Vdx, 


and  also  the  positive  pulse 


.+  * 


I.  -  /  (P.-Sdt  when  P, -1>0. 


'8) 


Here  t#  is  the  time  of  the  shock  wave  arrival  at  a  given  point; 

is  the  final  calculated  time.  The  dimensional  pulse  I  is  connected 

P 

with  the  dimensionless  pulse  Jp  by  the  relationship 


By  a  similar  method  it  is  possible  to  calculate  the  dimensionless 
pulses  of  the  velocity  pressure  Jy  by  the  formula 


d  X 

■ 


(9) 


where  the  connection  with  the  corresponding  dimensional  puise  will  be 


‘k  p  v2 


I*  *  /  “Y-  dt  -Jpt%  a*. 
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It  is  natural  that  the  integrals  in  expressions  (8)  and  (9)  should 
be  calculated  by  suitable  quadrature  formulas;  for  example,  during 
the  variable  step  of  integration  by  the  trapezoidal  rule.  In  this 
calculation  it  is  possible  to  also  determine  the  action  time  x+  for 
the  positive  pressure  phase,  using  the  equality 


-  x  -t. 


where  t  is  the  time  for  transition  of  the  excess  pressure  to  the 
negative  phase . 

During  the  numerical  solution  of  the  explosion  problem,  the 
pulses  in  the  cylindrical  and  spherical  cases  for  y  ■  1.4  were 
calculated  as  an  example .  The  variation  in  the  pulses  J*  and  Jy  and 
also  in  time  t+  and  7  that  depend  on  coordinate  R  for  v  ■  2  is 
presented  in  the  form  of  curves  in  work  [6].  Also  given  there,  for 
the  cylindrical  case,  are  curves  of  dimensionless  pressure  P,  which 
depends  on  t,  for  a  series  of  fixed  points  of  space. 

C.  Determining  the  Conneotion  Between 
Lagrange  and  the  Euler  Coordinates 

Let  us  designate  by  rQ  the  beginning  coordinate  for  a  particle 
of  gas.  Let  us  introduce  a  dimensionless  Lagrangian  coordinate  in 
the  form  of  n  *  (rQ/r°)v,  where  r°  Is  the  characteristic  length, 
determined  by  expression  (5)*  Then  the  connection  between  Lagrangian 
coordinate  n  and  the  Eulerian  coordinate  £  is  realized  with  the  help 
of  a  continuity  equation  in  this  fprm: 


By  using  the  tabulated  data,  it  is  possible  to  find  the  connection 
between  n  and  5  for  a  fixed  t  . 
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M.  USE  OP  THE  TABLES  IN  SOLVING  CERTAIN 
GAS-DYNAMIC  PROBLEMS  i 


A.  Determining  the  Parameters  of  a 
Point  Explosion  in  a>  Steady  Forward 
(fas  Flow 


Let  us  assume  that  an  explosion  occurs  in  a  gas  which  is  moving 
with  constant  velocity  U.  If  we  wish  to  examine  the  process  of  the' 
development  of  the  explosion  in  a  fixed  coordinate  system,  it  is 
necessary  to  recalculate  the  solution  by  using  the  Galileo-Newton 
transform.  Let  us  select  a  certain  fixed  system  of  Cartesian 
rectangular  coordinates  5^,  J£2»  *3  with  the  origin  in  the  center  of 
the  explosion  at  t  «  0.  Let  us  designate  by  U^,  U2>  projections 

'of  the  veldcity  vector  U  to  the  coordinate  axes,  and  by  0^,  ‘  02,  0^ 
the  corresponding  projections  of  the  velocity  vector  of  the  gas. 

1  1 

■  i 

Then  during  examination  of  the  gas  flow  in  the  fixed  coordinate 
system  Its  parameters  can  be  computed  by  using  the  conversion 
formulas 


*/  xj  ♦  V*  ‘V  "  vj  *  Ui 5 


here  Vj ,  p,  p  are  the  values  of  the  components  of  velocity,  pressure  , 
apd  density  in  the  coordinate  system  Xj  associated  with  the  forward 
moving  gas.  The  solution  of  the  problem  in  this  moving  coordinate 
system  corresponds  to  the  data  contained  in  the  tables.  Thus,  to 

determine  the  parameters  of  an  explosion  in  a  forward  gas  flow,  one 

♦ 

should  only  conduct  a  simple  recalculation  of  these  tabulated  data. 

When  conducting  such  a  recalculation  in  the  case  of  a  cylindrical 
or  a  spherical  explosion,  it  is  necfessary  to  take  into  account  the 
geometric  relationships  resulting  from  the  connection  between  Cartesian 
and  cylindrical  or  Cartesian  and  spherical  cdordinates.  ' 
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B.  Initial  Stage  of-  the  Reflection  of  a  Flat 
Blast  from  a  Parallel  Plane  Wall  and  of 
Cylindrical  and  Spherical  Waves  from  a 
Concentric  Cylindrical  or  Spherical 
Wall,  Respectively 

Let  us  examine  the  question  of  determining  the  pressure,  density 
and  velocity  of  a  gas  behind  the  front  of  a  reflected  shock  wave  for 
moments  of  time  close  to  the  moment  of  approach  of  the  wave  to 
a  wall  (designated  by  t#).  On  the  strength  of  the  assumed  geometry 
of  the  wall,  the  reflection  of  the  shock  wave  from  it  will  be  a  normal 
reflection,  but  the  flow  behind  the  front  of  the  reflected  wave  will 
be  one-dimensional  with  plane,  cylindrical,  or  spherical  symmetry. 

In  this  case,  the  pressure  p#,  the  density  p„  and  the  temperature  T„ 
in  the  reflected  wave  are  found  by  using  the  well-known  formulas  in 
gas  dynamics  (see,  for  example,  [7]): 


(10) 


L 

p 


JLA 


p  * 

(y  - 1 )  +i 

p_ 


L 

"p„  \ 


(id 


For  the  value  of  the  initial  velocity  of  the  reflected  wave  c„, 
we  have  the  expression 


■  r-  =  \  (12) 

i 

Since  the  wall  is  stationary,  on  the  strength  of  the  boundary  condition 
on  it  the  velocity  behind  the  front  of  the  reflected  shock  wave  w#  will 
be  small  for  moments  of  time  close  to  t*. 
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If  we  assign  the  distance  r#  from  the  point  of  impact  to  the 
wall  and  the  value  of  characteristic  length  r° ,  then,  knowing  the 
dimensionless  coordinate  R#  let  us  find,  according  to  the  value 
Rn  “  with  the  help  of  the  tables,  the  dimensionless  time  of  the 
approach  of  the  wave  to  the  wall  x#  and  by  it  time  t#,  and  also  the 
corresponding  values  of  q,  p  /p  and  p_/p_.  Furthermore,  by  using 
formulas  (10)-(12),  the  basic  parameters  of  the  reflected  shock  wave 
for  moments  of  time  close  to  t* . 

Let  us  note  that  formulas  ( 10 )  —  ( 12 )  for  the  problem  in  question, 
strictly  speaking,  are  accurate  only  at  the  moment  the  shock  wave 
directly  reflects  from  the  wall,  since  they  were  derived  for  a  uniform 
flow  behind  the  incident  wave.  However,  these  relationships  will 
be  approximately  fulfilled  for  those  moments  of  time  when  the 
variability  of  the  parameters  of  the  flow  behind  the  incident  wave  can 
be  disregarded. 

C .  Initial  Stage  of  Regular  Refleotion 
for  a  Plane ,  Cylindrical  or  Spherioal 
Blast  from  a  Flat  Surface 

Let  us  first  examine  the  plane  case.  Suppose  that  a  plane  blast 
strikes  an  absolutely  rigid  flat  wall  so  that  the  angle  between  the 
plane  of  the  shock  wave  and  the  wall  is  different  from  zero  and  equal 
to  a.  Then,  for  the  moments  of  time  close  to  the  moment  the  wave 
collides  with  the  wall,  from  relationships  on  the  shod'  wave  and  the 
boundary  conditions  on  the  wall  for  the  velocity  of  the  gas,  it  is 
possible  to  derive  the  analytical  dependences  between  the  parameters 
of  the  incident  and  reflected  waves  (see  [7]). 

Let  us  designate  by  o# ,  pt,  p§  respectively  the  reflection  angle, 
the  pressure  and  the  density  behind  the  reflected  wave  at  the  moment 
of  time  directly  after  the  reflection,  and  let  us  introduce  the 
following  parameters: 
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If  we  assign  angle  o  and  the  pressure  ratio  hw  in  the  moment  the  wave 
meets  the  reflecting  plane,  then  to  determine  we  have  the  quadratic 
equation 


(13) 


where 


Y-l  (l-kjtt 

I*  »  — — —  ;  M  c  - —  . 

Y+1  U|tk  ♦(it+k  )«* 

M  M 

The  value  of  the  relative  pressure  ht  behind  the  reflected  wave  is 
expressed  thusly: 


.  (in 

W,  -  *  (  |k  ♦  ID*  ) 

The  ratio  of  the  densities  pM/p„  and  of  the  temperatures  T-/T 

n  n 

can  be  found  with  the  help  of  formulas  (11).  Since  the  point  of 
intersection  of  the  incident  and  reflected  waves  moves  along  the 
plane,  the  velocity  of  the  reflected  shock  wave  c#  at  the  point  of 
reflection  can  be  obtained  from  the  relationship 


c  *  c  — - - •  \  / 

*  HR  01 

By  using  the  found  values  c*,  P*/p„  and  the  law  of  conservation 
of  mass  during  the  transition  through  the  explosion  surface,  it  is 
possible  to  determine  the  normal  velocity  component  of  the  gas  behind 
the  reflected  wave,  u#.  Viscous  component  is  found  by  the  condition 
of  its  continuum  during  the  transition  through  the  discontinuity. 

Let  us  note  that  the  lack  of  real  roots  of  equation  (13)  gives 
angular  values  to  which  the  irregular  (Mach)  reflection  sets  in. 
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Let  us  now  examine  the  case  of  the  reflection  of  a  cylindrical  or 
a  spherical  wave.  Suppose  that  an  explosion  takes  place  at  a  distance 
L  from  the  reflecting  plane  n  (in  the  cylindrical  case  for  simplicity 
let  us  consider  that  the  line  of  explosion  is  parallel  to  plajie  n). 

Let  us  designate  as  a  the  angle  at  which  a  cylindrical  or  a  spherical 
wave  will  approach  the  plane  II  at  a  certain  instant  t#  after  an 
explosion.  Prom  the  geometry  (Pig.  3),  it  is  evident  that  angle  o  is 
equal  to  the  angular  between  the  perpendicular  to  plane  n  and  the 
radius-vector  drawn  from  the  center  of  the  explosion  0  to  the 
reflection  point  of  the  wave  O'  at  the  instant  in  question,  t* 


0 


Pig.  3. 


By  applying  the  results  of  solving  the  problem  of  the  reflection 
of  a  plane  wave  to  the  reflection  of  an  element  of  a  cylindrical 
or  a  spherical  wave,  it  is  possible  to  determine  the  gas  parameters 
directly  at  the  moment  the  wave  element  reflects  from  the  plane. 
Having  designated  as  H  the  dimensionless  quantity  L/r°,  we  will  have 
the  equality 


CSlaa 


(16) 


Let  us  draw  our  interest  to  the  parameters  of  the  reflected 
wave  at  point  O'  of  plane  n  at  distance  rf  from  the  center  of  the 
explosion.  If  the  parameter  r°  is  known,  then  by  taking  Rt  *  r#/r° , 
it  is  possible  to  calculate  angle  a  by  formula  (16).  With  the  help 

of  the  tables  in  terms  of  the  value  R  ■  R«,  let  us  find  values 

n  * 
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q#,  t#  and  the  ratios  p^p,*,  If  with  the  derived  Pn/P00  angle 

a  corresponds  to  the  condition  of  a  regular  reflection,  then,  by 
using  the  designation  introduced  above  and  using  formulas  (11)  and 
(13)-(15),  let  us  calculate  the  parameters  of  the  reflected  wave  at 
instants  close  to  t# 

Furthermore,  suppose  that  the  dependence  between  Pn/Poo  and  the 
critical  angle  of  a  regular  reflection  is  known.  On  the  basis  of 
theoretical  analysis,  such  data  follow  from  formulas  (13)-(l6);  they 
are  given,  for  example,  in  book  [7].  Then,  for  a  given  value  L,  it 
is  possible  to  find  the  boundary  of  the  zone  a  regular  reflection  on 
plane  n. 

D.  Determining  the  Parameters  of  a  Flow 
During  an  Exploeion  on  a  Flat  Boundary 
Between  a  Gae  and  a  Solid  or  a 
Liquid  Medium 

Let  us  assume  that  the  upper  half-space  is  occupied  by  a  gas, 
and  the  lower  —  by  some  other  medium.  Let  us  examine  first  the 
idealized  case  when  the  lower  half-space  is  occupied  by  an  absolutely 
rigid  body.  To  determine,  with  the  help  of  the  tables  in  this  work, 
the  parameters  of  a  flow  during  the  explosion  of  a  charge  at  the 
boundary  of  the  gas  with  an  absolutely  solid  plane,  one  should  take 
the  value  |e°  in  all  cases  instead  of  the  value  E° . 

Suppose  that  the  lower  half-space  is  now  occupied  by  a  deformable 
medium,  but  the  flow  has  such  a  character  that  the  processes  occurring 
in  the  lower  half-space  very  weakly  affect  the  gas  flow  in  the  upper 
half-space.  Furthermore,  let  us  assume  that  the  time  of  the  shock 
wave  arrival  at  some  fixed  point  of  the  upper  half-space  is  known 
(for  instance,  this  time  is  measured  in  the  experiment).  As  such  a 
point  it  is  convenient  to  take  any  point  arranged  on  a  line, 
perpendicular  to  the  boundary  plane  and  passing  through  the  l.®  ter 
of  the  explosion  in  the  cylindrical  and  spherical  cases. 


( 

Then  It  is  possible  to  approximately  calculate  those  parts  of  the 
explosion  energy  E°  which  passed  into  the  upper  and  lower  half-spaces 
and  to  find  the  gas  flow  parameters  in  the  upper  half-space.  Actually, 
in  terms  of  the  coordinate  of  fixed  point  r#  and  the  time  of  the  shock 
wave  arrival  t*,  by  keeping  in  mind  formulas  (4),  (5)  and  by  using 
from  the  tables  the  dependence  of  R  upon  t«,  it  is  possible  to 
determine  the  energy  E®  which  was  released  into  the  upper  half-space 
It  is  evident  that  the  remaining  part  of  energy  E®  will  go  into  the 
lower  half-space,  since 

Bj°+  B*  -  E°.  (17) 

Knowing  energy  E®,  from  the  tables  it  is  possible  to  obtain  all  the 
gas-dynamic  parameters  of  interest  and  to  approximately  describe  the 
explosion  process  in  the  upper  half-space. 

Let  us  note  that  a  similar,  but  purely  experimental  approach 
to  finding  the  energy  distribution  between  the  two  half-spaces  was 
examined  in  work  [8]. 

E.  Plane  Explosion  at  the  Interfaoe  of 
Two  Identical  Gases  Having  Equal 
Initial  Pressures ,  but  Different 
Initial  Densities 

Let  us  discuss  the  question  of  using  the  tables  to  solve  the 
problem  of  determining  gas-dynamic  parameters  during  a  plane  explosion 
at  the  interface  of  the  gases  with  the  initial  parameters  y ,  p^, 
and  y,  pm,  p^. 

Let  us  assume  that  the  Lagranglan  coordinate  of  the  contact 
surface  is  n*  *  0.  Prom  the  condition  of  the  equality  of  the  pressures 
on  this  surface,  we  have 


PjW,  t)  -  Pt(0,  T)  , 


(18) 
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where  P1  and  P2  are  dimensionless  pressures  in  the  first  and  second  \ 
half-spaces,  respectively.  Let  us  take  now  as  functions  of  P1  and 
P2  the  dependences  obtained  for  a  homogeneous  medium.  Then  equality 
(18)  will  be  satisfied,  if  t1  ■  x2,  for  any  t.  Hence,  takinginto 
account  formulas  (5),  let  us  obtain  the  relation  between  the  fractions 
of  energy  that  went  out  into  first  half-space,  and  fractions  of 
E°  that  were  released  in  the  second  half-space,  and  namely: 


in  this  case,  formula  (17)  naturally  takes  place. 

Thus,  if  at  the  interface  of  identical  gases  there  occurred  a 
plane  explosion  with  energy  E°,  then  determining  the  physical 
characteristics  of  the  flow  in  both  half-spaces  can  be  done  from  the 
given  tables  while  taking  into  account  relationships  (17)-(19). 

There  is  a  possibility  of  broadening  the  circle  of  problems 
which  are  resolved  with  the  help  of  the  results  of  calculating  a 
point  explosion  in  a  homogeneous  gas.  Thus,  by  using  the  principle 
of  flat  cross  sections,  in  a  number  of  cases  it  is  possible  to 
approximately  find  the  parameters  of  an  explosion  in  a  heterogeneous 
medium.  If,  for  example,  in  the  case  of  a  cylindrical  explosion 
the  initial  density  is  slightly  changed  in  the  direction  of  the 
axis  of  the  explosion,  then  here,  on  the  strength  of  the  principle 
of  flat  cross  sections,  it  is  possible  to  conduct  calculations  for 
the  parameters  of  the  explosion  by  a  local  one-dimen3ional  theory. 
The  principle  of  flat  cross  sections  also  makes  it  possible  to  study 
the  explosion  along  a  line,  which  i$  different  from  a  straight  line, 
but  whose  curvature  is  everywhere  small.  Another  widely  known  and 
important  example  of  using  of  the  principle  of  flat  cross  sections 
will  be  examined  in  the  following  section. 


F.  Supplement  to  Problems  of  Hypersonic 
Flow  Around  Thin  Blunted  Plates  and 
Cylinders . 

Between  the  nonstationary  problem  of  an  explosion  and  the  problem 
of  steady  flow  over  bodies  at  hypersonic  velocities  a  known  analogy 
exists,  which  is  based  upon  the  principle  of  flat  cross  sections  [9]. 
Thus,  with  the  help  of  the  calculated  solution  for  plane  and 
cylindrical  explosions,  it  is  possible  to  approximately  determine  the 
hypersonic  flow  around  thin  plates  and  cylinders  with  small  blunting. 

In  this  case,  the  parameters  of  the  nonstationary  problem  are  replaced 
by  the  stationary  parameters,  with  t  and  E°  respectively  replaced  by 
xU^  and  ^xP0„uf  (xr/4  )v-1dv .  Here  the  following  designations  are  assumed 
x  —  the  coordinate  along  the  axis  of  the  body,  which  i3  measured  from 
its  tip;  y  -  the  coordinate  in  a  direction  perpendicular  to  the  surface 
of  the  plate  or  to  the  axis  of  the  streamlined  cylinder;  d  —  the 
characteristic  linear,  cross  section  size  of  blunt  part;  U  —  the 
velocity  of  the  incoming  gas  flow,  which  is  directed  along  the  axis  x; 
cx  ~  the  resist^1106  coefficient  of  blunt  part,  which  is  calculated 
per  unit  area  of  the  cross ■ section  of  the  body  and  referred  to  as  the 
velocity  pressure. 

During  such  a  change,  the  connection  between  the  dimensionless 
values  x/d,  y/d  and  x,  R  is  given  by  the  relationships 


where  ■  U^/a^  —  the  Mach  number  of  the  incident  flow.  Using 
formulas  in  (20),  it  is  possible  with  the  help  of  the  given  tables 
to  determine  the  gas  flow  parameters  near  a  plate  and  a  cylinder  at 
sufficient  distance  from  the  blunt  portion.  Specifically,  it  is 
possible  to  obtain  the  pressure  distribution  on  body  p/p^,  the  shape 
of  the  leading  shock  wave  (i.e.,  the  dependence  of  y/d  on  x/d)  and 
the  values  of  the  basic  gas-dynamic  functions  behind  it. 
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The  relationships  given  above  in  sections  B  and  C  (see  also  [7]) 
make  it  possible  to  calculate  the  parameters  for  the  gas  directly 
behind  the  front  of  the  reflected  shock  wave  which  appears  upon  the 
regular  reflection  of  the  leading  shock  wave  (being  formed  in  front 
of  the  body)  from  a  flat  rigid  surface.  Such  a  reflection  can  take 
place  during  a  hypersonic  flow  over  a  body  In  channels  or  during  the 
motion  of  thin  blunt  bodies  in  a  gas  near  solid  boundaries. 
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